Current response of ac-driven nanoelectromechanical systems in single-electron tunneling regime 
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We investigate electric current in a single-electron tunnelling device weakly coupled to an ac-driven under- 
damped harmonic nanomechanical oscillator. In the linear regime, the current can respond to the external fre- 
quency in a resonant as well as in an anti-resonant fashion. The main resonance is accompanied by an additional 
resonance at a half of the external frequency. 

PACS numbers: 85.85.+j,73.23.Hk,62.40.+i 



I. INTRODUCTION 

During last decade the nanoelectromechanical systems 
(NEMS) were widely studied because of richness of phys- 
ical effects and their broad spectrum of the practical 
functionalities!. With small inertial masses and high fre- 
quencies of the nanomechanical resonators, together with the 
ultra-sensitive mechanical displacement detection capabilities 
of electronic devices used as detectors, NEMS have great 
promise for application in the force microscopy and mass- 
sensing 26 . Many of the applications currently envisioned, 
such as switches, relays, and actuators, are directly related to 
the effect of mechanical motion on the electric properties of 
the system. 

Recently, the first observation of quantum superposition of 
states of a mechanical resonator has been reported 7 . Bottom- 
up fabricated NEMS devices are prospective candidates for 
investigation of quantum effects in NEMS, since they have 
frequencies in GHz range, and these frequencies are tun- 
able by the gate voltage 8,9 . These devices, made of sus- 
pended carbon nanotubes, function in single-electron tunnel- 
ing (SET) regime. The coupling between the charge of the 
SET and mechanical degrees of freedom is provided by the 
position dependence of the capacitance between suspended 
beam and the underlying gate, and of the tunnel rates between 
the reservoirs and the SET island. It was theoretically pre- 
dicted previously 10 that for certain energy-dependent tunnel 
rates even at weak coupling between the SET and mechanical 
motion the latter can considerably influence electric current 
— the phenomenon of strong mechanical feedback. At strong 
coupling (see Ref. 11), the effect of the mechanical motion is 
even stronger, and the physics is dominated by polaronic ef- 
fect. This leads to multistabilities and switching in the regime 
of a>o <K F, where wo and T are the resonator frequency and 
the typical tunnel rate, respectively-^. For wo » F, one ob- 
serves Franck-Condon effec t 13,14 . 

Recently, careful experiments on electric transport in sus- 
pended carbon nanotubes have been performed 15 . In addition 
to source, drain and gate electrodes, an antenna emitting ac 
radiation was introduced close to the vibrating beam. A vari- 
ety of interesting results have been obtained, including reso- 
nances, switching events, and gate voltage dependence of the 
quality factor of the mechanical resonator. There are indica- 
tions that strong coupling regime may have been observed, at 
least in some parameter range. 



In this Article, we consider simultaneous effects of ac and 
stochastic forces on transport in SET-based NEMS. A sim- 
ilar problem was previously considered in Ref. 16. Our 
manuscript is different from Ref. [l6j in several respects: We 
use a different theoretical formulation. In contrast to Ref. [Is, 
we consider the energy dependent tunneling rates (in addition 
to the energy dependence given by Fermi functions). We also 
calculate the dependence of the electric current on the driv- 
ing frequency. We find that the regime of strong mechanical 
feedback, which was discussed in reference 10, disappears al- 
ready at weak ac force. The main effect of the ac force is to 
create resonances when the external frequency coincides with 
wo and wo/2. These resonances can enhance as well as sup- 
press (anti-resonance) the current. We find that in some cases 
a resonance at o>o can be accompanied by an anti-resonance at 
wo/2 and vice versa. 



II. MODEL AND THEORY 

The system we study is an SET device coupled to a har- 
monic oscillator with the frequency a>o> mass M and the qual- 
ity factor Q » 1, The most obvious specific realization is a 
suspended beam (carbon nanotube) 8 , see Fig. Q] The beam is 
suspended between two leads; the bias voltage V;, - Vl-Vr 
is applied. The beam thus at the same time is an island where 
the electron charge is quantized. 

The island is capacitively coupled to an underlying gate at 
which voltage V g is applied. This coupling originates from the 
dependence of the equilibrium position of the oscillator on the 
charge state of the island and is characterized by the force F s , 
that acts on the oscillator from the gate electrode when there 
is an extra electron on it. The coupling can be quantified by 
a dimensionless constant A - F^/hMoj^, which is the relative 
shift of the oscillator energy resulting from a single tunnelling 
event in respect to zero point motion energy. We assume that 
oscillator in addition is driven by an external antenna which 
creates an ac signal with the amplitude Vd and the frequency 
a). Fig. | l(b)| shows the equivalent circuit: The dot represents 
the island which is connected to four voltage sources via the 
capacitances. We assume that the displacement of the island x 
is much smaller then the distance between the island and the 
gate electrode d\ and both these distances are smaller then the 
distance between the island and the antenna di {x <K d\ <K 
dT). It follows that the gate capacitance is much bigger than 
the capacitance between the island and the antenna Co •^ C g . 
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FIG. 1: (a) Picture of the setup: a beam is suspended over two electrodes can oscillate freely. The beam is capacitively coupled to the gate 
electrode. Bias voltage V b = V L — Vr is applied to the electrodes, and the beam is driven by the ac signal applied to the antenna, (b) Schematic 
version of setup: the dot in the center corresponds to the beam which is connected to four voltage sources via corresponding capacitances. 



The electrostatic energy of this system in the first order of 
displacement can be written in the following way, 
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where E ch = e 2 /(C L + C R + C° g + C D ) and q = C L V L + C R V R + 
CgVg + CdVd are the charging energy and the induced charge 
on the island by the electrodes respectively (C° is gate capac- 
itance for x — 0). n represents the charge state of the island, 
and F(n) is the force acting on the island from all four elec- 
trodes, including the driving force. Charge transfer process 
at a given n is characterized by an energy difference between 
final and initial state, which is the difference of electrostatic 
energy given by Eq. ([TJ plus (minus) an energy cost associ- 
ated with addition (extraction) an electron to (from) the corre- 
sponding electrode, 

A£r" +1 = E ch in + i - @) - eV L - F st (n)x; 
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where F st = F(n + 1) - F(ri). 

In what follows, we generalize the treatment of Ref. UM 
to the case when the oscillator is ac driven. At weak cou- 
pling, the motion of the oscillator is classical, and our qualita- 
tive evaluation of the electron transport in this regime is based 
on the master equation for the distribution function P„(x, v, t), 
where v is the velocity of the oscillator. We assume that 
the bias is low enough so that only two charge states, n and 
n + 1, are important for the transport. This distribution func- 
tion obeys the Boltzmann equation 17 , 
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Here T is the total force acting on the oscillator which is the 
sum of the elastic force, the friction force with the damping 
coefficient y and the electric forces acting from the electrodes. 
The last ones include the external driving force which we 
choose to be periodic in time. The "collision integral" St [P] 
consists of two terms: the first representing the SET and the 
second describing interaction of the island with the environ- 
ment with temperature T em , There are four tunnel rates, T^ s , 
where the subscripts L and R denote tunneling through the left 
or right junction, and the superscripts + and - correspond to 
the tunneling to and from the island, respectively. In Eq. ((5]) 
F* = F^ +F^. Each rate is a function of the corresponding ad- 
dition energy A£^ R associated with the addition/removal of 
an electron to/from the island. 

We can simplify Eq. ([3]) if we restrict ourselves to the fast 
tunnelling regime a> ~ «o ^ E,, where E, = E + + E~ is to- 
tal tunnelling rate, which means that the oscillator coordinate 
varies so slowly that E(x) hardly changes between two succes- 
sive tunneling events. In this regime we can apply the adia- 
batic approximation, arriving to the Kramers equation 10 , 
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where T\ is the total force T excluding the friction force. Co- 
efficients A(x), D(x) and k{x) are the following functions of 
the tunneling rates, 
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Here D(x) is the diffusion coefficient which contains the term 
describing the driving of the oscillator by the stochastic force 
Fst, In further calculation we assume that environment tem- 
perature sets the lowest energy scale in our system, kbT em , <K 
min(£ , c / 1 , eV/,), so that we drop first term in Eq. (0. Further- 
more, k(x) is the drift coefficient and if this is positive the cor- 
responding term can be interpreted as a dissipation. For some 
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FIG. 2: (a) The I-V characteristics of the system for the fixed value of the bias voltage V b = 0.35mV and the amplitude of the driving signal 
Vrf = 0.3 V. The blue (red) curve corresponds to the system out of (at) the resonance. Inset shows the stability diagram, the current in main 
figure is calculated along the blue line and two red points indicate the parameter values, which were used for calculating the resonances in 
(b). (b) The current response on the driving frequency for the different values of the quality factor. The upper panel shows the response in 
the current enhancement region V g = -0.031V (resonance) and the lower shows the response in the current suppression region V g = -0.03 V 
(unti-resonance) 



tunneling rates and sufficiently high quality factors the drift 
coefficient can be negative and in this case it gives accumula- 
tion of energy. The third term A(x) has a rather trivial effect, 
it can be included into the oscillator potential energy and just 
renormalizes the elastic force, not leading to qualitative new 
effects. Thus we disregard this term in further consideration. 

In the single electron transport regime between n and n + 1 
the island is subject to the force which is the sum of two terms: 
the total electric force F(n) and the stochastic force F st due to 
the stochastic electron transfers events. The solution of the 
equation of motion for the oscillator subject to such force can 
be written in the following way: x = Xo + 5x, where xq is 
due to the total electric force and is a non-stochastic function 
obtained from the solution of the equation of motion. The 
second term is a stochastic variable which originates from the 
stochastic force. If we assume that the quality factor Q is 
sufficiently big (2^1) an d we are out of the resonance (wo + 
ai), then the last term can be parameterized by the energy E 
and the phase as follows, ~\/2E/Mcl>() sin <p. 

The master equation can be solved for two limiting cases. 
The first is when xq » 6x, dependence of the diffusion and 
feedback coefficients on 6x can be disregarded. After averag- 
ing Eq. (O over the phase we obtain the following, 
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to the second regime which corresponds to xq <s 6x. This 
second regime has been studied in Ref. 10. The stationary 
distribution function in this regime is given by 



P 2 (E) = P' exp 
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Here D(E) - (D(x) cos 2 <f>)^, and k{E) - (k(x) cos 2 <J))q, where 
the angular brackets denote the average over the oscillation 
period, {A(y)) - J(d(p/27T)A( V2esin0). For the entire energy 
range, whereas we can not solve the Kramers equation exactly, 
we can approximate the solution as follows, 



P(E) = 



Pi(E) if E < E Q ; 
Pi(E) if E > E . 



(13) 



The average electric current flowing through the system can 
be calculated by using this distribution function 
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We can derive the time averaged distribution function which 

is 



Pi(£) = P exp 



(k(Xq)), . 

'(D(xo)}, 1 
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Here (k(xq)}, and {D(xq)}, are the time averaged drift and dif- 
fusion coefficients. If for some values of parameters the drift 
coefficient becomes negative, the probability increases which 
means that the oscillator will be driven to higher amplitudes 
until the energy E hits the value Eq, where the system switches 



Using these equations it is possible to calculate the current 
flowing through the system if we know the energy dependence 
of the tunneling rates. 



III. NUMERICAL RESULTS 

For illustration we have chosen the system considered in 
Ref. [l5l which is similar to the setup depicted in Fig. Q] 
The exponential energy dependence of the tunneling barriers 
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FIG. 3: The dependence of the current on the driving frequency for the fixed values of voltages. Each panel correspond to three different values 
of gate values: (a) - V g = -0.0315V, (b) - V g = -0.0296V, (c) -V g = -0.03V. Red points on the insets show parameter values, which were 
used to calculate the resonances in the main figures. 



is typical for such SET system 18 , so we have considered the 
following rates, 



F + LR = 2b L , R e- a ^U M AEt R ) ; 
r t ,« - b L , R e- a ^ E l«f F (AE L , R ) , 



(16) 



where a^ R and b^ R are some constants, and /f(E) - 1/[1 + 
exp(E/ki,T e „ v )] is the Fermi function. The factor 2 accounts 
for the spin degeneracy of the state n. For more concrete- 
ness, we take C L = C R = 160aF and C g = 1000C D = 80aF. 
The temperature of the transported electrons we choose to be 
T env = 200mK and the value of the oscillator frequency a>o we 
take to be 1 GHz. 

Fig. |2(a)| represents the current-voltage characteristics for 
the above described system in the single electron transport 
regime. For the fixed value of bias voltage V/, = 0.35mV we 
have varied the gate voltage and have calculated the current. 
The blue curve on Fig. |2(a)| corresponds to the immovable is- 
land. When the island is in resonance with the driving signal, 
which we take with the amplitude Vd = 0.3V, the current is 
modified due to the oscillations (red curve). For the different 
gate voltages the current at the resonance is either suppressed 
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FIG. 4: The current- voltage characteristics of the system for two dif- 
ferent values of temperature of tunneling electrons: red correspondes 
to T = 0.2K and blue-7 = 1A*. Inset showes dependence of current 
on the temperature of the electrons for fixed values of the electrod 
voltages: V g = -0.0315V, V b = 0.35mV and V d = 0.3V 



or enhanced. The linear response of the current on the driving 
frequency is shown on the figure [2(b)| for the different values 
of the quality factor. Depending on which region we take, the 
current enhancement or the current suppression, we get the 
linear resonance or anti-resonance peaks. 

The electric force F(n) in Eq. (Q]), acting on the island from 
the electrodes can be decomposed as the sum of three terms: 

F(n) = MV, N)cos 2 (cjt) + f 2 (V,N)cos(cjt) + f 3 (V,N) (17) 

where fi(V,N), fi(V,N) and f$(V,N) depend on the charge 
state of the island and the voltages applied on the electrodes 
(bias, gate and driving voltages). Because of the first term in 
Eq. ( fTTb the amplitude of the oscillations increases at the half 
of the driving frequency a> - a>o/2. This gives an additional 
response peak of the current at this value of the frequency. 

Fig. [3] shows the current as the function of the driving 
frequency for the fixed values of voltages {V/, - 0.35mV, 
Vd = IV). Three different panels correspond to the different 
values of the gate voltage. Thus, by changing the gate voltage 
it is possible to get either two resonance peaks (c), two anti 
resonance peaks (b) or one resonance and one anti-resonance 
peaks (a) at the resonance frequency and at it's half. 

Increasing the temperature of transferred electrons causes 
broadening of coulomb picks and coulomb diamonds. Fig. @] 
shows current dependence on gate voltage for different values 
of electron temperature. For higher values of the tempera- 
ture the current pick decreases in value and broadens (blue 
curve). Inset represents temperature dependence of the cur- 
rent for fixed values of applied voltages, which shows that for 
high temperatures the height of Coulomb peak decreases ex- 
ponentially. 



IV. CONCLUSION 

In this Article, we studied the behavior of the "self- 
detecting" NEMS in the regime were it can modeled as a sin- 
gle electron transistor weakly coupled to a harmonic mechan- 
ical oscillator. In particular, we studied the dependence of 
electric current on the external driving frequency. We demon- 
strated that current is modified when the system is at the res- 
onance. We also found that since current can be enhanced 



by mechanical motion as well as suppressed, depending on 
the applied voltages, the driving of the resonator can result in 
the resonance peak (associated with the current enhancement) 
as well as in an antiresonance dip (associated with the cur- 
rent suppression). Additionally, the driving force is propor- 
tional contains doubly-oscillating terms, which results in the 
appearance of an additional resonance peak/dip in the spec- 
trum at the half of the natural frequency coq. Depending on 
the applied voltages, both resonances can be peaks, dips, or 
one can have a peak and a dip. The latter situation is realized 
provided the current is enhanced by mechanical motion at low 
amplitudes and suppressed at high amplitudes, or vice versa. 
The regions of strong backaction discussed in Ref. 10 for the 
system studied above appear to be for high, where the system 
is not in the single electron tunnelling regime anymore, and 
our simple analysis can not be applied to this situation. 

This Article treats the simplest situation — when the oscil- 
lator is linear, and the coupling is weak. Experiments 15 clearly 
indicate that in certain voltage and power range the resonator 



can be driven into non-linear and/or strong coupling regime 19 . 
With the increase of the driving power the oscillator enters 
non-linear regime, and bistabilities and switching effects be- 
tween states with different amplitude have been experimen- 
tally observed. Current response to the external frequency in 
a non-linear system will have a more complicated form that 
discussed in this Article, in particular, resonances at the fre- 
quencies multiple of wo/2 are expected. These effects require 
a more extensive analysis which will be performed elsewhere. 
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